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INTRODUCTION 

The metal-insulator transition (MIT) in doped semiconductors with a random 
distribution of impurities is a unique quantum phase transition in the sense that both 
disorder and electron-electron interaction play a key role (see for example Refs. 1 and 2). 
The metallic phase of the transition is characterized by a finite electrical conductivity 
at T = 0, while the conductivity in the insulating phase vanishes in the limit of zero 
temperature. From a theoretical point of view, the correlation length in the metallic 
phase and the localization length in the insulating phase diverge at the critical point 
with the same exponent u, i.e., they are proportional to \N/N C — l\~ u in the critical 
regime of the MIT, and the value of v provides important information about the MIT. 
Here, N is the impurity concentration and N c is the critical concentration for the MIT. 
Since direct experimental determination of v is extremely difficult, researchers have 
usually determined, instead of u, the value of the conductivity critical exponent /i 
defined by 

a(0) = a*(N/N c -ir (1) 

immediately above N c (0 < N/N c — 1 <C 1). Here, cr(0) is the conductivity extrapolated 
to T = and a* is a prefactor. Values of v are then obtained assuming the relation 
v — jj, [3] which is valid for three-dimensional systems without electron-electron in- 
teraction. With a number of nominally uncompensated semiconductors fi ~ 0.5 has 
been obtained [2]. One of the best studied nominally uncompensated semiconductors is 
Si:P. Rosenbaum et al. studied both the doping- induced MIT and the uniaxial-stress- 
induced MIT, and showed that Eq. (1) describes <j(0) of Si:P with a single exponent 
fi pa 0.5 over a wide range 0.001 < N/N c — 1 < 1 [4]. Here, uniaxial stress was used for 
"tuning" because fine control of N/N c — 1 is difficult. In the same material, however, 
fi pa 1.3 was claimed about ten years later for a narrow regime N c < N < 1.1N C by 



a different group [5]. Moreoever, fj, pa 1 was reported recently for the MIT driven by 
uniaxial stress [6]. As the origin of the discrepancy, inhomogeneous distribution of the 
impurities has been pointed out [7]. 

For the case of melt- (or metallurgically) doped samples, which have been employed 
in most of the previous studies including Refs. 4 — 6, the spatial fluctuation of N due 
to dopant striations and segregation can easily be on the order of 1% across a typical 
sample for the four-point resistance measurement that has a length of ~5 mm or larger 
(see for example Ref. 8). For this reason, it will not be meaningful to discuss physical 
properties in the critical regime (e.g., \N/N C — 1| < 0.01), unless one evaluates the 
macroscopic inhomogeneity in the samples and its influence on the results. In order to 
rule out the ambiguity arising from the inhomogeneity, we prepared 70 Ge:Ga samples by 
neutron-transmutation doping (NTD) of isotopically enriched 70 Ge single crystals. The 
NTD method inherently guarantees the random distribution of the dopants down to 
the atomic level [9,10]. We show from the conductivity measurements at T = 0.02 — 1 K 
that tr(O) of the NTD 70 Ge:Ga samples is described by Eq. (1) with fj, = 0.50 ± 0.04 
over a wide range 4 x 10~ 4 < N/N c - 1 < 0.4 [11]. 

In order to determine v without assuming /j, — u, we have analyzed the tem- 
perature dependence of the conductivity on the insulating side of the MIT in the 
context of variable-range-hopping (VRH) conduction within the Coulomb gap [12]. 
Low-temperature (T < 0.5 K) conductivity of the insulating 70 Ge:Ga samples obeys 
<t(T) = cr(T) exp[— (To/T) 1 / 2 ], which is predicted by the VRH theory [13], with an 
appropriate temperature dependence in the prefactor <r(T). Magnetic field and tem- 
perature dependence of the conductivity of the samples are subsequently measured in 
order to determine directly the localization length £ and the impurity dielectric suscep- 
tibility Ximp as a function of N in the context of the theory. This kind of determination 
of £ and Ximp was performed for compensated Ge:As by Ionov et al. [14] They found 
£ oc (1 - N/N c )~ v and Ximp oc (N c /N - with v = 0.60 ± 0.04 and C = 1-38 ± 0.07, 
respectively, for samples having N up to 0.96iV c . The significance of their result is 
the experimental verification of the relation 2is pa ( that had been predicted by scal- 
ing theories [15]. However, the critical exponents of compensated samples are known 
to be different from those of nominally uncompensated samples [11]. Therefore, our 
determination of £ and xi mp in nominally uncompensated samples is important. 

The previous effort to measure the impurity dielectric susceptibility as a function 
of N has also contributed greatly. Hess et al. found ( = 1.15 ± 0.15 in nominally 
uncompensated Si:P [16]. Since fj, ~ 0.5 was determined for the same series of Si:P 
samples [4], the relation 2/x ~ ( was again valid. Katsumoto has found ( ~ 2 and /i ~ 1 
for compensated Al . 3 Ga .7As:Si, i.e., again, 2/i pa ( applies [17]. Thus, in these cases 
the conclusion 2is pa ( was reached indirectly, by assuming n — v. We, on the other 
hand, determine v directly, i.e., we do not have to rely on the assumption /j, — v in 
order to study the behavior of the localization length £ near the MIT. 

According to theories [2] on the MIT which take into account both disorder and 
electron-electron interaction, the critical exponents do not depend on the details of 
the system, but depend only on the universality class to which the system belongs. 
Moreover, there is an inequality v > 2/3, which is expected to apply generally to disor- 
dered systems irrespective of the presence of electron-electron interaction [18]. Hence, 
if one assumes \i = v, which is derived for systems without electron-electron interaction, 
/x pa 0.5 violates the inequality. This discrepancy has been known as the conductivity 
critical exponent puzzle. Kirkpatrick and Belitz have claimed that there are logarithmic 
corrections to scaling in universality classes with time-reversal symmetry, i.e., when the 
external magnetic field is zero, and that \i pa 0.5, found at B = 0, should be inter- 
preted as an "effective" exponent which is different from a "real" exponent satisfying 
/i > 2/3 [19]. Therefore, comparison of ji with and without the time-reversal symme- 



try, i.e., with and without external magnetic fields becomes important. We study the 
MIT of 70 Ge:Ga in magnetic fields up to B = 8 T and show that \x changes from 0.5 
at B = to 1.1 at B > 4 T [20]. The same exponent // = 1.1 is also found in the 
magnetic-field-tuned MIT for three different samples, i.e., 

a(N,B,T^0)cx[B c (N)-Bf, (2) 

where B C (N) is the critical magnetic field for concentration N. Moreover, an excellent 
finite-temperature scaling [2] 

a(N, B, T) oc T x f(\B c (N) - B\/V), (3) 

where x/y is equivalent to //, is obtained with the same value of //. The phase diagram 
on the (N, B) plane is successfully constructed, and we find a simple scaling rule which 
a(N,B,T — > 0) would obey and derive from a simple mathematical argument that 
/i = jj! as has been observed in our experiment. 



EXPERIMENTS 
Sample preparation 

All of the Ge:Ga samples were prepared by neutron-transmutation doping (NTD) 
of isotopically enriched 70 Ge single crystals. The basic idea of NTD is as follows. 
Suppose that a nucleus in a crystal of a semiconductor captures a thermal neutron. 
After the capture, the nucleus is not necessarily stable. If it is stable, the element 
remains unchanged, but if it is not, it decays and transmutes into a new element which 
may act as a dopant. This is NTD. Practically, a crystal is placed in a nuclear reactor 
which produces thermal neutrons. Since the neutron field produced by a reactor is large 
enough to guarantee a homogeneous flux over the crystal dimensions and the small 
capture cross section (typically 10~ 24 cm 2 ) of semiconductors for neutrons minimizes 
"self-shadowing" , NTD is known to produce the most homogeneous, perfectly random 
distribution of dopant down to the atomic level [9]. 

As for Ge, there are five stable isotopes: 70 Ge (20.5%), 72 Ge (27.4%), 73 Ge (7.8%), 
74 Ge (36.5%), and 76 Ge (7.8%). The numbers in the parentheses represent the natural 
abundance. These five stable isotopes of Ge undergo the following nuclear reactions 
after capturing a thermal neutron. 

Ga (EC, 11.2 d), (4) 

(5) 
(6) 

As (/3, 82.8 min) , (7) 
As (/?, 11.3 hrs) -> USe ({3, 38.8 hrs) . (8) 

Here, Jn is a neutron, EC and (5 denote electron capture and f3 decay, respectively, and 
the time in the parentheses represents the half life. Note that the natural Ge form 
both acceptors (Ga) and donors (As and Se) after NTD. Empirically, the impurity 
compensation is known to affect the value of the critical exponent [11]. To avoid the 
impurity compensation, we use isotopically enriched 70 Ge. 

The Czochralski grown, chemically very pure 70 Ge crystal has isotopic composition 
[ 70 Ge]=96.2 at. % and [ 72 Ge]=3.8 at. % [10]. The as-grown crystal is free of dislocations, 
p type with a net electrically-active-impurity concentration less than 5x 10 11 cm -3 . The 
thermal neutron irradiation was performed with the thermal to fast neutron ratio of 
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Figure 1. Electrical conductivity as a function of T 1 / 2 for NTD 70 Ge:Ga. From bottom to top in 
units of 10 17 cm- 3 , the Ga concentrations are 1.853, 1.856, 1.858, 1.861, 1.863, 1.912, 1.933, 2.004, 
2.076, 2.210, 2.219, 2.232, 2.290, 2.362, and 2.434, respectively. 

«30:1. The small fraction of 72 Ge becomes Ge which is stable, i.e., no further acceptors 
or donors are introduced. The post-NTD rapid thermal annealing at 650 °C for 10 sec 
removed most of the irradiation-induced defects from the samples. The short annealing 
time is important in order to avoid the redistribution and/or clustering of the uniformly 
dispersed 71 Ga acceptors. The concentration of the electrically active radiation defects 
measured with deep level transient spectrometry (DLTS) after the annealing is less 
than 0.1% of the Ga concetration, i.e., the compensation ratio of the samples is less 
than 0.001. The dimension of most samples used for conductivity measurements was 
6x0.9x0.7 mm 3 . Four strips of boron-ion-implanted regions on a 6x0.9 mm 2 face 
of each sample were coated with 200 nm Pd and 400 nm Au pads using a sputtering 
technique. Annealing at 300 °C for one hour activated the implanted boron and removed 
the stress in the metal films. 

The concentration of Ga acceptors after NTD is determined from the time of 
thermal-neutron irradiation. The concentration is proportional to the irradiation time 
as long as the same irradiation site and the same power of a nuclear reactor are em- 
ployed. 



Low-temperature measurements 



The electrical conductivity measurements were carried out down to temperatures 
of 20 mK using a 3 He- 4 He dilution refrigerator. All the electrical leads were low-pass 
filtered at the top of the cryostat. The sample was fixed in the mixing chamber and a 
ruthenium oxide thermometer [Scientific Instrument (SI), RO600A, 1.4x1.3x0.5 mm 3 ] 
was placed close to the sample. To measure the resistance of the thermometer, we used 
an ac resistance bridge (RV-Elekroniikka, AVS-47). The thermometer was calibrated 
against 2Ce(N0 3 )3-3Mg(N0 3 )2-24H 2 (CMN) susceptibility and against the resistance 
of a canned ruthenium oxide thermometer (SI, RO600A2) which was calibrated com- 
mercially over a temperature range from 50 mK to 20 K. We employed an ac method 
at 21.0 Hz to measure the resistance of the sample. The power dissipation was kept 
below 10~ 14 W, which is small enough to avoid overheating of the samples. The output 
voltage of the sample was detected by a lock-in amplifier (EG&G Princeton Applied 
Research, 124A). All the analog instruments as well as the cryostat were placed inside a 
shielded room. The output of the instruments was detected by digital voltmeters placed 
outside the shielded room. All the electrical leads into the shielded room were low-pass 
filtered. The output of the voltmeters was read by a personal computer via GP-IB 
interface connected through an optical fiber. Magnetic fields up to 8 T were applied in 
the direction perpendicular to the current flow by means of a superconducting solenoid. 



RESULTS AND DISCUSSION 

Temperature dependence of conductivity in metallic samples and the critical 
exponent for the zero-temperature conductivity 

The temperature dependence of the electrical conductivity mostly for the metallic 
samples is shown in Fig. 1. The temperature variation of the conductivity of disordered 
metal is governed mainly by electron-electron interaction at low temperatures [1], and 
can be written as 

Aa(T) = a{T) - <r(0) = mVf , (9) 

where 

Here, F is a dimensionless and temperature-independent parameter characterizing the 
Hartree interaction and D is the diffusion constant [1], which is related to the conduc- 
tivity via the Einstein relation 

a = (dn/dn)e 2 D, (11) 

where (dn/dfj,) is the density of states at the Fermi level. In various reports such 
as Refs. 5 and 10, er(0) was obtained by extrapolating cr(T) to T = assuming y/T 
dependence based on Eq. (9). One should note, however, that it is sound only in 
the limit of Acr(T) <C er(0) ~ cr(T), where D can be considered as a constant, i.e., 
m is constant, and that the inequality is no longer valid as N approaches N c from 
the metallic side since <r(0) also approaches zero. In such cases m in Eq. (9) is not 
temperature independent and Acx(T) may exhibit a temperature dependence different 
from y/T. To examine this point in our experimental results, we go back to Fig. 1. 
We see there that A<r(T) of the bottom five curves are not proportional to y/T : while 
Aa(T) of the other higher N samples are described by cx y/T. The close-ups of cr(T) for 
the six samples with positive da/dT in the scale of y/T and T 1 / 3 are shown in Figs. 2(a) 
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Figure 2. Conductivity as a function of (a) T 1 / 2 and (b) T 1 / 3 , respectively, near the metal-insulator 
transition. From bottom to top in units of 10 17 cm -3 , the concentrations are 1.853, 1.856, 1.858, 1.861, 
1.863, and 1.912, respectively. The upper and lower dotted lines in each figure represent the best fit 
using the data between 0.05 K and 0.5 K for the first and the third curves from the top, respectively. 
Each fit is shifted downward slightly for easier comparison. 



and 2(b), respectively. The upper and lower dotted lines represent the best fit using 
the data between 0.05 K and 0.5 K for the samples with N = 1.912 x 10 17 cm" 3 and 
N = 1.861 x 10 17 cm" 3 , respectively. Each fit is shifted downward slightly for the sake 
of clarity. From this comparison, it is clear that a T 1 / 3 dependence rather than a \JT 
dependence holds for samples in the very vicinity of the MIT. The opposite is true for 
the curve at the top. This means that the \/T dependence in Eq. (9) is replaced by a 
yi/3 dependence as the MIT is approached. 

A T 1 / 3 dependence close to the critical point for the MIT was predicted originally 
by Al'tshuler and Aronov [21]. They considered an interacting electron system with 
paramagnetic impurities, for which they obtained a single parameter scaling equation. 
At finite temperatures, they assumed a scaling form for conductivity according to the 
scaling hypothesis: 

*=j^fiWLr), (12) 



where £' is the correlation length and Lt = yhD/ksT is the thermal diffusion length. 
When L T > fi(£'/L T ) = A + A 1 £'/L T , which is equivalent to Eq. (9). In the 
critical region, where L T <C £' — > oo, Eq. (12) should be reduced to 



e 2 



" = C W < 13 > 

Combining this equation and Eq. (11), they obtained a oc T 1 / 3 . The T 1 / 3 dependence 
was also predicted from numerical calculations that consider solely the effect of disor- 
der [22]. 

It is not clear whether the argument in Refs. 21 or 22 is applicable to the present 
system or not, but there is an experimental fact that VT dependence of conductivity 
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Figure 3. (a) Conductivity a as a function of (T/er) 1 / 2 . From bottom to top in units of 10 17 cm" 3 , the 
concentrations are 1.858, 1.861, 1.863, and 1.912, respectively. The solid lines denote the extrapolation 
for finding tr(0). (b) Zero-temperature conductivity cr(0) vs the dimensionless distance N/N c — 1 from 
the critical point on a double logarithmic scale. The dotted line represents the best power-law fit by 
cr(0) oc (N/N c - iy where fi = 0.50 ± 0.04. 
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Figure 4. The logarithm of the conductivity as a function of T -1 / 2 for insulating samples, the 
concentrations from bottom to top in units of 10 17 cm" 3 arc 1.717, 1.752, 1.779, 1.796, 1.805, 1.823, 
1.840, 1.842, 1.843, 1.848, 1.850, 1.853, 1.856, and 1.858, respectively. 



changes to T 1 / 3 as N c is approached. It is important that we find a consistent method 
that allows the determination of cr(0) for both the cases. For this purpose, we follow 
Al'tshuler and Aronov's manipulation [21] of eliminating m and D in Eqs. (9)-(ll) and 
obtain 

a(T) = a(0)+m'yjT/a(T), (14) 

where w! = eA^J (dn/dp) , which is temperature independent. In the limit of Aer(T) <C 
<j(0) ~ u{T) , this equation gives the same value of <x(0) as Eq. (9) does. When 
cr(0) <C °~(T), it yields a T 1 / 3 dependence for a(T). Thus, it is applicable to both 
a/T and T 1 / 3 dependent conductivity. From today's theoretical understanding of the 
problem, Eqs. (9) and (14) are valid only for L T 3> and their applicability to the 
critical region is not clear, because the higher-order terms of the (3 function [23] which 
were once erroneously believed to be zero do not vanish [24]. Nevertheless, we expect 
Eq. (14) to be a good expression for describing the temperature dependence of all 
metallic samples because it expresses both \f f and T 1 / 3 dependences as limiting forms. 
Then, based on Eq. (14), we plot o~(T) vs ^T/a(T) for the four close to N c samples 
in Fig. 3(a). The data points align on straight lines, which supports the adequacy of 
Eq. (14). The zero-temperature conductivity <r(0) is obtained by extrapolating to T = 0. 
The curve on the top of Fig. 3(a) is for the sample with the lowest N among the ones 
showing \/T dependence at low temperatures, i.e., this sample has the largest value of 
A<t(T)/(t(0) among \/T samples. The value of <r(0) obtained for this particular sample 
using Eq. (14) differs only by 0.6% from the value determined by the conventional 
extrapolation assuming Eq. (9). This small difference is comparable to the variation 
arising from the choice of the temperature range in which the fitting is performed. 
Therefore, the extrapolation method proposed here is compatible with the conventional 
method based on the oc \/T extrapolation. 

Based on this analysis the MIT is found to occur between the first and second 
samples from the bottom in Fig. 3(a), i.e., 1.858 x 10 17 cm~ 3 < N c < 1.861 x 10 17 cm~ 3 . 
Thus, unlike the case for Si:P [5,7], N c is fixed already within an accuracy of 0.16% 
and the evaluation of the critical exponent /i will not be affected by the ambiguity in 
the determination of N c . Figure 3(b) shows the zero-temperature conductivity <r(0) as 
a function of N/N c — 1 on a double logarithmic scale. A fit of Eq. (1) (dotted line) 
is excellent all the way down to (N/N c — 1) = 4 x 10 -4 . The fitting parameters are 
H = 0.50 ± 0.04, N c = (1.860 ± 0.002) x 10 17 cm" 3 , and a* = (40 ± 2) S/cm. We note 
that pi = 0.46 ± 0.18 is obtained even when we use only the four samples closest to N c 
for the fitting. 

Variable-range-hopping conduction in insulating samples and the critical 
exponents for localization length and dielectric susceptibility 

The temperature dependence of the conductivity <r(T) of insulating samples is 
shown in Fig. 4. The electrical conduction of doped semiconductors on the insulating 
side of the MIT is often dominated by variable-range hopping (VRH) at low tempera- 
tures. The temperature dependence of a(T) for VRH is written in the form of 

a(T)=a(T)exp[-(T /Tn (15) 

where p — 1/2 for the excitation within a parabolic-shaped energy gap (Coulomb gap), 
and p = 1/4 for a constant single-particle density of states around the Fermi level [13]. 
The temperature dependence of <t(T) contributes greatly to the temperature depen- 
dence of cr(T) near N c because the factor T /T in the exponential term becomes very 
small, i.e., the temperature dependencies of the prefactor and that of the exponential 
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Figure 5. Conductivity multiplied by T 1 / 3 vs (a) T x l 2 and (b) T 1 / 4 . From bottom to top in 
units of 10 17 cm~ 3 , the concentrations are 1.848, 1.850, 1.853, and 1.856, respectively. 



term become comparable. Theoretically, cr(T) is expected to depend on temperature as 



a(T) oc T 



(16) 



but the value of r including the sign has not been derived yet for VRH conduction with 
both p — 1/2 and p = 1/4. 

As we have seen in Fig. 2(b), the temperature variation of the low-temperature 
conductivity of the 70 Ge:Ga samples within ±0.3% of N c is proportional to T 1 / 3 . Since 
both the T 1 / 3 dependence of the conductivity and the VRH with p = 1/2 are results of 
the electron-electron interaction in disordered systems, they can be expressed, in princi- 
ple, in a unified form. Moreover, the electronic transport in barely metallic samples and 
that in barely insulating samples should be essentially the same at high temperatures 
so long as the inelastic scattering length and the thermal diffusion length are smaller 
than, or at most comparable to the correlation length or the localization length. So, the 
temperature dependence of conductivity at high temperatures should be the same on 
both sides of the transition. Such behavior is confirmed experimentally in the present 
system, i.e., as seen in Fig. 2(b) the conductivity of samples very close to N c shows a 
T 1 / 3 dependence at T « 0.5 K, irrespective of the phase (metal or insulator) to which 
they belong at T = 0. Based on this consideration we fix r = 1/3. Figure 5 shows oT~ r 
with r = 1/3 for four samples (N/N c = 0.993, 0.994, 0.996, and 0.998) as a function 
of (a) T -1 / 2 and (b) T -1 / 4 . All the data points lie on straight lines with p — 1/2 in 
Fig. 5(a) while they curve downward with p = 1/4 in Fig. 5(b). This dependence is 
maintained even when we change the values of r between r = 1/2 and 1/4. Thus we 
conclude that the conductivity of all samples on the insulating side for N up to 0.998^ 
is described by the theory for the VRH conduction where the excitation occurs within 
the Coulomb gap, i.e., Eq. (15) with p = 1/2. 

Based on these findings, we evaluate the N dependence of T in Eq. (15) with 
p — 1/2 and r = 1/3. Figure 6 shows To as a function of 1 — N/N c . The vertical 
and horizontal error bars have been estimated based on the values of T obtained with 
r = 1/2 and r = 1/4, and the values of 1 - JV/(1.858 x 10 17 cm~ 3 ) and 1 - JV/(1.861 x 
10 17 cm -3 ), where 1.858 x 10 17 cm -3 is the highest concentration in the insulating phase 
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Figure 6. T determined by a(T) oc T 1 / 3 exp[— (To/7 1 ) 1 / 2 ] as a function of the dimensionless 
concentration 1 — N/N c . 

and 1.861 x 10 17 cm" 3 is the lowest in the metallic phase, respectively. According to 
theory [13], T in Eq. (15) is given by 

in SI units, where e(N) is the dielectric constant, and £(iV) is the localization length. 
Here, we should note that the condition T < To is needed for the VRH theory to be 
valid, i.e., T has to be evaluated only from the data obtained at temperatures low 
enough to satisfy the above condition. This requirement is fulfilled in Fig. 6 for all the 
samples except for the one with N = 0.998iV c . Concerning this latter sample, we will 
include it for the determination of 7 (Fig. 8) and £ and Ximp (Fig. 9) but not for the 
calculation of the critical exponents. 

Our next step is to separate T into e and £ in the framework of the theory of 
VRH conduction with and without a weak magnetic field [13]. For £/A <C 1, the 
magnetoconductance is expressed as 
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where A = yh/eB is the magnetic length in SI units. According to Eq. (18), the 
magnetic-field variation of In a at T = const, is proportional to B 2 , i.e., 



- In a(B, T) = — In <x(0, T) + C 2 (T)B\ 



(19) 



and the slope C 2 (T) in the above equation is proportional to T~ 3 / 2 . In order to demon- 
strate that these relations hold for our samples, we show for the iV = 0.989iV c sam- 
ple — ]na(B,T) vs B 2 in Fig. 7(a) and C^T) determined by least-square fitting of 
-diner /dB 2 vs T~ 3/2 in Fig. 7(b). Since Eq. (18) is equivalent to 



7 = C 2 (T)/T- 3 / 2 w 0.0015 (e/h) 2 ^T 3/2 , 



(20) 



£ is given by 
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Figure 7. (a) Logarithm of cr _1 vs B 2 at constant temperatures for the sample having N = 1.840 x 
10 17 cm~ 3 . From top to bottom the temperatures are 0.095 K, 0.135 K, and 0.215 K, respectively. 
The solid lines represent the best fits, (b) Slope —d\na/dB 2 vs T~ 3 / 2 for the same sample. The solid 
line represents the best fit. 
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Figure 8. Coefficient 7 defined by Eq. (20) as a function of T . 




N c /N-1 




Figure 9. (a) Dielectric susceptibility Ximp arising from the impurities vs N c /N— 1. (b) Localization 
length f vs 1 - N/N c . 



In this way we have determined 7 as a function of T for nine samples (Fig. 8). The 
value of 7 is almost independent of T , and if one assumes the form of 7 oc Tj/, one 
obtains a small value of 5 = 0.094 ± 0.005 from least-square fitting. We determine £ 
and Xi m p = e — €h from Eqs. (17) and (21), and show them in Fig. 9 as a function of 
1 — N/N c and N c /N — 1, respectively. Here, e h is the dielectric constant of the host 
Ge, and hence, Ximp is the dielectric susceptibility of the Ga acceptors. We should note 
that both £ and Ximp are sufficiently larger than the Bohr radius (8 nm for Ge) and 
€h = 15.4 [25], respectively. According to the theories of the MIT, both £ and Ximp 
diverge at N c as f (JV) oc (1 - N/N c )- V and x imp (iV) oc {N c /N respectively. We 

find, however, both £ and Ximp do not show such simple dependencies on N in the range 
shown in Fig. 9, and that there is a sharp change of both dependencies at N as 0.99iV c . 
On both sides of the change in slope, the concentration dependence of £ and x imp are 
expressed well by the scaling formula as shown in Fig. 9. Theoretically, the quantities 
should show the critical behavior when N is very close to N c . So v = 1.2 ± 0.3 and 
( = 2.3 ± 0.6 may be concluded from the data in 0.99 < N/N c . However, the other 
region (0.9 < N/N c < 0.99), where we obtain v = 0.33 ± 0.03 and C = 0.62 ± 0.05, is 
also very close to N c in a conventional experimental sense. 

As a possible origin for the change in slope, we refer to the effect of compensa- 
tion. Although our samples are nominally uncompensated, doping compensation of less 
than 0.1% may be present due to residual isotopes that become n-type impurities after 
NTD. In addition to the doping compensation, the effect known as "self compensa- 
tion" may play an important role near N c [26]. It is empirically known that the doping 
compensation affects the value of the critical exponents [11]. Rentzsch et al. studied 
VRH conduction of n-type NTD Ge in the concentration range of 0.2 < N/N c < 0.91, 
and showed that T vanishes as T oc (1 — N/N c ) a with a ~ 3 for K = 38% and 
54%, where K is the compensation ratio [27]. Since a ~ v + ( [Eq. (17)], we find for 
our NTD 70 Ge:Ga samples a = 3.5 ± 0.8 for 0.99 < N/N c < 1 and a = 0.95 ± 0.08 
for 0.9 < N/N c < 0.99. Interestingly, a = 3.5 ± 0.8 agrees with a m 3 found for 
compensated samples. Moreover, we have recently proposed the possibility that the 
conductivity critical exponent fi ~ 1 in the same 70 Ge:Ga only within the very vicinity 
of N c (up to about +0.1% of N c ) [28]. An exponent of fi = 0.50 ± 0.04, on the other 
hand, holds for a wider region of iV up to 1.4iV c as we have seen in Fig. 3(b). Again, 
/i ~ 1 near N c may be viewed as the effect of compensation. Therefore, it may be 



^ I I I ■■MIT— ~~ I I I I I 

0.0 0.2 0.4 0.6 0.8 




Figure 10. Conductivity of the sample having N — 2.004 x 10 17 cm 3 as a function of T 1 / 2 at 
several magnetic fields. The values of the magnetic induction from top to bottom in units of tesla are 

0. 0, 1.0, 2.0, 3.0, 4.0, 4.7, 5.0, 5.3, 5.6, 6.0, 7.0, and 8.0, respectively. 

possible that the region of iV around N c where v « 1 and /i ~ 1 changes its width 
as a function of the doping compensation. In the limit of zero compensation, the part 
which is characterized by v ~ 1 and /j, ~ 1 vanishes, i.e., we propose v = 0.33 ± 0.03, 
C = 0.62 ± 0.05, and fi = 0.50 ± 0.04 for truly uncompensated systems and that the 
relation /i — v [3] is not satisfied. In compensated systems, on the other hand, \i = v 
may hold as it does in the very vicinity of N c . However, the preceding discussion needs 
to be proven experimentally in the future by using samples whose compensation ratios 
are controlled precisely and systematically. 

Metal-insulator transition in magnetic fields 

Figure 10 shows the temperature dependence of the conductivity of the sample hav- 
ing N = 2.004 x 10 17 cm -3 for several values of the magnetic induction B. Application of 
the magnetic field decreases the conductivity and eventually drives the sample into the 
insulating phase. This property can be understood in terms of the shrinkage of the wave 
function due to the magnetic field. In strong magnetic fields and at low temperatures, 

1. e., when g/ie-B ^> ksT, the conductivity shows another \/T dependence [1] 




(22) 




Figure 11. Zero-temperature conductivity a(N,B,0) vs normalized concentration n = 

[a(N, 0, 0)/er*(0)] 2 '° = N/N c (0) - 1, where er(7V,0,0) is the zero-temperature conductivity and a*(0) 
is the prefactor both at B = 0. From top to bottom the magnetic induction increases from 1 T to 
8 T in steps of 1 T. The dashed curve at the top is for B = 0. The solid curves represent fits of 
a(N, B, 0) oc [n/n c {B) - l]^ B \ For B > 6 T, we assume \i = 1.15. 




where 

mB = Ti^t/r If -i^iir- < 23 > 

One should note that Eqs. (9) and (22) are valid only in the limits of [a(N,0,T) — 
a(N, 0, 0)] < a(N, 0, 0) or [a(N, B, T) - a(N, B, 0)] < a(N, B, 0). It is for this reason 
that we have observed a T 1 / 3 dependence rather than the \/T dependence at B — 
in Fig. 2 as the critical point [cr(N, 0, 0) = 0] is approached from the metallic side. 
However, Fig. 10 shows that the \/T dependence holds when B ^ even around the 
critical point. Hence, we use Eq. (22) to evaluate the zero-temperature conductivity 
a(N,B,0) in magnetic fields. Since m B is independent of B, the conductivity for 
various values of B plotted against y/T should appear as a group of parallel lines. This 
is approximately the case as seen in Fig. 10 at low temperatures (e.g., T < 0.25 K). 

The zero-temperature conductivity a(N,B,0) in various magnetic fields obtained 
by extrapolation of a(N,B,T) to T = based on Eq. (22) is shown in Fig. 11. Here, 
a(N, B,0) is plotted as a function of the normalized concentration: 

n^[a(iV,0,0)K(0)] 2 - . (24) 

Since the relation between N and cr(N, 0, 0) was established in Fig. 3(b) as <j(N, 0, 0) = 
a*(0)[N/N c (0) - l] - 50 where JV c (0) = 1.860 x 10 17 cm" 3 and <r*(0) = 40 S/cm, n is 
equivalent to iV/iV^O) — I. Henceforth, we will use n instead of N because employing 
n reduces the scattering of the data caused by several experimental uncertainties, and 
it further helps us concentrate on observing how cr(N, B, 0) varies as B is increased. 
Similar evaluations of the concentration have been used by various groups. In their 
approach, the ratio of the resistance at 4.2 K to that at a room temperature is used 
to determine the concentration [5]. The dashed curve in Fig. 11 is for B = 0, which 
merely expresses Eq. (24), and the solid curves represent fits of 

a(N,B,0) = a (B)[n/n c (B) - l]« B l (25) 



The exponent n{B) increases from 0.5 with increasing B and reaches a value close to 
unity at B > 4 T. For example, // = 1.03 ± 0.03 at B = 4 T and fi = 1.09 ± 0.05 at 




Figure 12. Zero-temperature conductivity a(N,B,0) vs magnetic induction B. From bottom to 
top, the normalized concentrations denned by Eq. (24) are 0.04, 0.09, and 0.22, respectively. 



B = 5T. When B > 6 T, three-parameter [a (B), n c (B), and fJ,(B)] fits no longer give 
reasonable results because the number of samples available for the fit decreases with 
increasing B. Hence, we give the solid curves for B > 6 T assuming fi(B) = 1.15. 

We show cr(N, B, 0) as a function of B in Fig. 12 for three different samples. When 
the magnetic field is weak, i.e., the correction Aa B (N, B, 0) = a(N,B,0) — a(N, 0, 0) 
due to B is small compared with a(N, 0, 0), the field dependence of Aub(N, B, 0) looks 
consistent with the prediction by the interaction theory [1], 

In larger magnetic fields, a(N, B, 0) deviates from Eq. (26) and eventually vanishes 
at some magnetic induction B c . For the samples in Fig. 12, we tuned the magnetic 
induction to the MIT in a resolution of 0.1 T. We fit an equation similar to Eq. (25), 

v(N,B,0) = a' (n)[l - B/B c {n)]^ n \ (27) 

to the data close to the critical point. As a result we obtain /x' = 1.1 ± 0.1 for all of the 
three samples. The value of // depends on the choice of the magnetic-field range to be 
used for the fitting, and this fact leads to the error of ±0.1 in the determination of //. 
In Fig. 13 we show that // = 1.1 yields an excellent finite-temperature scaling [Eq. (3)]. 
Note that the data both on the metallic side and on the insulating side are included in 
this scaling plot. Here we employ B c obtained by fitting Eq. (27), x — 1/2 from the 
fact that \J~T dependence holds in magnetic fields even around the critical point, and 
y — x/fj,' — 0.45, i.e., none of them are treated as a fitting parameter. Hence, Fig. 13 
strongly supports fj/ — 1.1. 

From the critical points n c (B) and B c (n), the phase diagram at T = is con- 
structed on the (N, B) plane as shown in Fig. 14. Here, n c (B) for B > 6 T shown by 
triangles are obtained by assuming \x = 1.15. The vertical solid lines associated with the 
triangles represent the range of values over which n c (B) have to exist, i.e., between the 
highest n in the insulating phase and the lowest n in the metallic phase. Solid diamonds 
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Figure 13. Finite-temperature scaling plot for the field-induced metal-insulator transition in the 
70 Ge:Ga sample having n — 0.09. 
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Figure 14. Phase diagram of 70 Ge:Ga at T = 0. The solid circles and the open triangles represent the 
critical concentrations n c , and the solid diamonds and the open boxes represent the critical magnetic 
induction B c . 



represent B c for the three samples in which we have studied the magnetic-field-induced 
MIT. Estimations of B c for the other samples are also shown by open boxes with error 
bars. The boundary between metallic phase and insulating phase is expressed by a 
power-law relation: 

n = C 3 B 13 . (28) 

From the eight data points denoted by the solid symbols, we obtain C 3 = (1.33 ±0.17) x 
KT 3 T-P and (3 = 2.45 ±0.09 as shown by the dotted curve. The shift of N c in magnetic 
fields was studied theoretically by Khmel'nitskii and Larkin [29]. They considered a 
noninteracting electron system starting from 

a(N,B,0)^^f 2 (B a 'a, (29) 

where £' is the correlation length. They claimed that the argument of the function f 2 
should be a power of the magnetic flux through a region with dimension This means 

a(N,B,0)^^f 2 (e/X), (30) 

where A = is the magnetic length, and hence, a' = 1/2. In order to discuss the 

shift of the MIT due to the magnetic field, they rewrote Eq. (30) as 

a(iV, J B,O)^|^0(tA 1 ^), (31) 

based on the relation in zero magnetic field 

e oc t~ v . (32) 

Here, t is a measure of distance from the critical point in zero field, e.g., 

t=[N/N c (0)-l}. (33) 

The zero point of the function gives the MIT, and the shift of the critical point for 
the MIT equals 

N C (B) - N c {0) oc B 1/2u . (34) 

Thus, (3 = l/(2z/) results. In the present system, however, this relation does not hold, 
as long as we assume /j, — v [3]. Experimentally, we find [3 = 2.5, while l/(2z/) = 
1/(2//) = 1 for 70 Ge:Ga at B = 0. 

Based on the phase diagram we shall consider the relationship between the two 
critical exponents: /x for the doping-induced MIT and // for the magnetic-field-induced 
MIT. Suppose that a sample with normalized concentration n has a zero-temperature 
conductivity a at B ^ and that \n/n c (B) - 1] < 1 or [1 - B/B c (n)} < I. From 
Eqs. (25) and (27), we have two expressions for a: 

a = a (n/n c - 1)" (35) 

and 

a = a' (l-B/B c f. (36) 
On the other hand, we have from Eq. (28) 

n/n c = (B/B c )-P = [1 - (1 - B/B c )}^ »! + /?(!- B/B c ) (37) 
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Figure 15. Normalized zero-temperature conductivity a(N,B,0)/a' (n) and <r( AT, B, 0)/[/3 M cto(S)] 
as functions of [1 — B/B c (n)] and [n/n c (2?) — l]//3, respectively, where /3 = 2.5 and ^ = 1.1. The solid 
line denotes a power-law behavior with the exponent of 1.1. The open and solid symbols represent the 
results of the magnetic-field-induced metal-insulator transition (MIT) in the range (1 — B/B c ) < 0.5 
for three different samples (n = 0.04, 0.09, and 0.22) and the doping-induced MIT in constant magnetic 
fields (4, 5, 6, 7, and 8 T), respectively. 



in the limit of (1 — B/B c ) <C I. This equation can be rewritten as 

(n/n c -l)//3 » (1-B/B C ). 
Using Eqs. (35), (36), and (38), we obtain 

</ (l - B/B c y' w /3"<7„(1 - B/B c y. 
Since Eq. (39) has to hold for arbitrary B, the following relations 

a' = P»ct 

and 



(38) 

(39) 

(40) 
(41) 



fi = fi 
are derived. 

In Fig. 15 we see how well Eq. (41) holds for the present system. We have already 
shown in Fig. 12 that // = 1.1 ± 0.1 is practically independent of n. Concerning the 



exponent //, however, its dependence on B has not been ruled out completely even 
for the highest B we used in the experiments. This is mainly because the number 
of available data points at large B is not sufficient for a precise determination of fi. 
In Fig. 15 the results of the doping- induced MIT for B > 4 T (solid symbols) and 
the magnetic-field-induced MIT for three different samples (open symbols) are plotted. 
Here, we plot a(N, B, 0)/[^ a (B)} vs [n/n c (B) - l]/(3 with (3 = 2.5 and /i = 1.1 for 
the doping-induced MIT, and a(N, B,0)/<Jo(B)' vs [1 — B/B c (n)} for the magnetic- 
field-induced MIT. Figure 15 shows that the data points align exceptionally well along 
a single line describing a single exponent // = // = 1.1. 

We saw in Fig. 11 that /i apparently takes smaller values in B < 3 T, which 
seemingly contradicts the above consideration. We can understand this as follows. We 
find that the critical exponent fi in zero magnetic field is 0.5 which is different from 
the values of \i in magnetic fields. Hence, one should note whether the system under 
consideration belongs to the "magnetic-field regime" or not. In systems where the MIT 
occurs, there are several characteristic length scales: the correlation length, the thermal 
diffusion length, the inelastic scattering length, the spin scattering length, the spin-orbit 
scattering length, etc. As for the magnetic field, it is characterized by the magnetic 
length A = ^Jh/eB. When A is smaller than the other length scales, the system is in 
the "magnetic-field regime." As the correlation length £' diverges at the MIT, A < £' 
holds near the critical point, no matter how weak the magnetic field is. When the 
field is not sufficiently large, the "magnetic- field regime" where we assume fi = 1.1 to 
hold, is restricted to a narrow region of concentration. Outside the region, the system 
crosses over to the "zero-field regime" where fi = 0.5 is expected. This is what is seen 
in Fig. 11. 

CONCLUSION 

We have measured the electrical conductivity of NTD 70 Ge:Ga to study the metal- 
insulator transition, ruling out an ambiguity due to inhomogeneous distribution of 
impurities. The critical exponent fi ~ 0.5 in zero magnetic field for doped semicon- 
ductors without impurity compensation has been confirmed. On the insulating side 
of the MIT, while the relation 2v xs ( predicted by scaling theories [15] holds for 
0.9 < N/N c < 1, the critical exponents for localization length and impurity dielectric 
susceptibility change at N/N c as 0.99. The small amount of doping compensation that 
is unavoidably present in our samples may be responsible for such a change in the ex- 
ponents. We have also measured the conductivity in magnetic fields up to B = 8 T 
in order to study the doping-induced MIT (in magnetic fields) and the magnetic-field- 
induced MIT. For both of the MIT, the critical exponent of the conductivity is 1.1, 
which is different from the value 0.5 at B — 0. The change of the critical exponent 
caused by the applied magnetic fields supports a picture in which \i varies depending 
on the universality class to which the system belongs. The phase diagram has been 
determined in magnetic fields for the 70 Ge:Ga system. 
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